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Abstract 

In this study, variations with respect to temperature of the incre- 
ments of enthalpy AH and Gibbs energy AG, arising in the dissolution 
of proteins in water, have been investigated by the methods of statis- 
tical thermodynamics. In this formalism, effective electric field E and 
total dipole moment M are taken as thermodynamical variable. In 
obtaining the free energy the partition function given by A. Bakk, 
J.S. Hoye and A. Hansen; Physica A, 304, (2002), 355-361 has been 
used in a modified form. In the constructed semi-phenomenological 
theory, the experimental data are taken from the study of Privalov 
^ [2| and the relevant parameters have been found by fitting to the 
experimental curves. The variations of the increments of enthalpy 
AH and Gibbs energy AG have been investigated in the temperature 
range 265 — 350i^. 
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1 Introduction 



Proteins is the common name of the complex macro molecules which are 
formed by gathering of a great number amino-acids and which play unremit- 
ting roles in the life times of all of the living creatures. There are approxi- 
mately twenty different amino-acids which take place in the structure of the 
proteins and which are chain molecules. The order or arrangement of this 
twenty units is sensitive and definite. This order, determines the nature of 
the protein and as a result identifies its function [3]. Each protein molecule 
is formed by on original combination and consecutive placing of different 
numbers of amino-acids. 

Proteins under their own natural physiological conditions, spontaneously 
open out from one dimensional structures to three dimensional structures in 
the time intervals starting from milliseconds and getting to minutes 13] . 
The three dimensional structures of the proteins depend on various environ- 
mental factors in the cell. Unfolded proteins return to their natural states 
when the environmental factors are removed exhibiting the feature that the 
three dimensional structure of a protein solely depends on its amino-acid or- 
dering information. This situation gives the hope to some research workers 
that protein folding could be reconstructed by computer simulation j^J. But 
the difficulty of this simulation even for small proteins is in view. 

In this study, the model belonging to the proteins is presented in section 2. 
In section 3, the thermodynamical quantities concerning the variation of the 
increments of additional enthalpy AH and Gibbs energy AG with respect 
to temperature have been expressed and these parameters are related to 
the parameters of the proteins with the help of statistical thermodynamics. 
At the end, from the curves belonging to the increments of enthalpy and 
Gibbs energy, the conclusions relevant to the structures of the proteins are 
presented. 

2 Model 

A simple model of protein unfolding has been developed by Bakk Bakk 
has started the model by choosing the water molecules as classical electric 
dipoles. he has demanded that non polar dissolution should include impor- 
tant physical properties in order to model the effect on the folding of the 
proteins. In the model used, in the medium of the formed electric field, the 
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folding and unfolding of the proteins have been taken as the basis. Here, the 
electric field is not an external field but has been used in the modelling of 
the ice-like behavior exhibited by the water molecules around the non polar 
surfaces. The electric field is a result of the effective behavior of the non 
polar dissolvent applied to the protein unfolding. 

The modified partition function of the protein system could be taken as 

m 

Z(T, E, = exp[exp(^)z(T, E, N = 1)]. (1) 

Here, \x represents the chemical potential and the term £ = exp(/5/z) 
corresponding to the pH of the system has been introduced to the partition 

function. z(T,E,N = 1) = 47rexp~^~2 — ^-sinh/3e e indicates single particle 
partition function. Where effective energy is given bye e = e + b<m> |7j . 
On the other hand average dipole moment has the form < m >= L[j3e e ). 

3 Variations of the Enthalpy and Gibbs En- 
ergy Increments with Respect to Temper- 
ature 

Taking electric field E and total dipole moment M as thermodynamical vari- 
ables, the first law of thermodynamics could be written for the proteins in 
the form: 

dU = TdS - MdE. (2) 
In this case, change in the enthalpy is; 

dH = TdS + EdM. (3) 

Enthalpy and entropy increments could also be written by taking the 
total differentials of enthalpy H = H(S, M) and entropy S = S(T, M) which 
leads to: 

dH =(^s)MdS+( M ) s dM (4) 

P) Q Q 

dS=(g f ) M dT+(—) T dM. (5) 
After substituting Eq.(4) and Eq.(5) in Eq.(3) 
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has been obtained. 

Writing down the expressions T = (§§)m, E = (f§)s and (|| ) M = 
in Eq.(6) and then integrating one could get for the enthalpy increment: 

rT r M ( T ) <99 

AH= C M dT+ { E + T(—) T ]dM. (7) 

JT t J M(T t ) OM 

Using the first law of thermodynamics for the proteins which has been 
given by Eq.(2), change in the Gibbs energy 

dG = -SdT + EdM (8) 

could be written down. With a similar approach, taking the total differentials 
of the Gibbs energy G = G(T, M) and total dipole moment M = M(T, E) 
and then substituting then in Eq.(8) together with the expressions E = 
(fjH)r and —S = (§^)m leads after integration to the expression which gives 
the Gibbs energy increment: 

,T ,M ( T) 

AG = - / SdT + / EdM. (9) 

JT t J M(T t ) 

The free energy of the system has been obtained from the relation F = 
—kTlnZ giving: 

Aexp(-(3(\bm 2 - ^))7rsinh(/3e e ) 

F = We • (10) 

The thermodynamical quantities of the system has been determined using 



the entropy relation S = — which in turn leads to: 



S = 2ex P(~^ 2 - ^ x { (_2^cosh(^ e ) + (11) 

P^e 

(4 + (3bm 2 - 20 n) sinh(/5e e ))}. 



The total dipole moment has been calculated from the expression M = 
(§§)t resulting in the following: 

M = n ' K * 2 ™- x {e e /?cosh(/3e e ) - sinh(/3e e )}. (12) 

P 
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The energy of the system is given by U = — \^ and when the partition 
function is substituted one obtains: 

u = - "» x { (-4« e ,aco S h(,a fe ) + 

2n{2 + (3bm 2 - 2/3/i) sinh(^e e ))}. (13) 



Additional heat capacity at effective field has been calculated by making 
use of the equation AC E = —k/3 2 (^) E which gives: 

exp(—/3(^bm 2 — a))kn 

(4e e /3(f3bm 2 + 2/3/i)) cosh(/3e e ) + (8 + /3(6 2 m 4 /3 - 8/i)) sinh(,5e e ) - 

f34bm 2 (pfi - 1) sinh(/3e e ) + 4f3 2 (e 2 e + /i 2 ) smh(/3e e )}. (14) 



Using the expression AC M - AC E = -T(§§ ) T [(p )m] 2 which has been 
obtained for additional heat capacity at constant total dipole moment; 

= exp{-P{\bm 2 - f,))kn ^ { _^ 2 + ^ _ ^ cosh(/9ee) + 

[8 + /3(& 2 m 4 /3 - 8/i - 46m 2 (/9 / u - 1) + 4/3(e 2 + ^ 2 ))] sinh(/3e e ) - 

k 2 [3 2 



x [e e( 5(4/36m 2 - cosh(/3e e ) - 



(-2/3e e cosh(/?e e ) + (2 + (3 2 e 2 e ) sinh(/3e e )) 

(4 + [3(bm 2 + 2{3e e - 2fi)) sinh(/3e e )] 2 }(15) 



has been calculated 

,r I 



Another quantity {§ti)t when calculated gives: 



O t = ~ k(3ee - (16) 

When Eq.(15) and Eq.(16) are substituted in Eq.(7), after integration gives 
for the enthalpy increment: 
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AH = ^Z^gj x {r , CTp( ^WH-5)+^ ) x 
[2e e (bm + t - ke + {k - l)e e ) cosh(^) + fc(2(fc - l)T t (e - e e ) - 

bm(2T 4 + me e ) + 2e e /i) sinh(^)] + Texp( ^ tH> ) x 

[-2e e (6m + e - ke + (k - l)e e ) cosh(-^-) + fc(-2(fc - l)T(e - e e ) + 

kT 

bm{2T + me e ) - 2e ef i) sinh(^)] }. (17) 

Similarly, substituting Eq.(ll) into Eq.(9), after integration and writing 
down Eq.(12) in this expression, one obtains for the Gibbs energy increment: 

1/—rW ftgggHgh fern! , 2m 

AG = P[ 2 2fcTtT J x {T t exp( T h Tt ) x 

e e 2« 

[-ee e cosh(^-) + fcT t (e + e e ) sinh^)] - 

fcm 2 I 2/i 

{Texp( ^i_X) x [-ee, cosh(^) + kT{e + e e ) sinh(^)] }. (18) 

Now the numerical applications of the enthalpy and Gibbs energy incre- 
ments will be taken under consideration. The unfolding and then folding of 
the protein molecules take place by a thermodynamical mechanism. In other 
words, the original structures of the macromolecules are determined ther- 
modynamically. In the following figures, dashed lines represent experimental 
study where as solid lines theoretical results. The physical quantities e, b and 
\i of the graphs which have been obtained in this semi phenomenological the- 
ory are determined by fitting to the experimental results of Privalov PJ E] . In 
Figs. (1) and (2) the variations of the thermodynamical functions enthalpy 
and Gibbs energy increments with respect to temperature are represented 
for two of the selected proteins. Yellow (belonging to Myoglobin) and red 
(belonging to Ribonuclease) solid lines indicate the variations with respect 
to temperature of increments of the thermodynamical quantities which have 
been obtained theoretically. 
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Figure 1: The variation of the additional enthalpy with temperature. 




Figure 2: The variation of the additional Gibbs energy with temperature. 

4 Conclusions 

The enthalpy increments AH between the folding and unfolding states, has 
a negative value in the cold folding region. This behavior shows that cold 
folding has been realized by giving heat out and this is negative heat capacity, 
on the other hand in the AH > case the system is in the warm folding 
region and the heat capacity is positive. 

The greatest contribution to the natural structure are provided by hy- 
drophobic interactions. In an unfolding chain structure, hydrophobic amino- 
acid side chains are unfolded by a regular water molecule lattice and the 
entropy of the system decreases. But when the chain undergoes folding and 
hydrophobic side chains are collected at intermolecular regions by getting a 
way from water, the regular water molecules will be free and as increase in 
the total entropy of the system takes place. 

In each protein molecule AG Gibbs potential increment, is a fraction of 
the quantity kT room in the folding and unfolding states at each degrees of 
freedom and AG Gibbs potential increment is not a monotonic function of 
temperature. Gibbs free energy increment between the folding and unfolding 



states has a maximum that is AG > and thus the unfolded form is stable 
and at both sides of this maximum this increment becomes negative. This 
means that two transitions take place. One of them is the unfolding of the 
proteins at high temperature where as the other is the folding of the proteins 
at low temperature. This phenomenon is known as cold denaturation in 
other words cold unfolding. As a result, in the cold denaturation transition 
hidden heat is negative and just opposite to this in the warm denaturation 
transition hidden heat has a positive value. 

During the folding of the protein entropy increases that is to say in re- 
turning to the original structure of the protein it is not preferred from the 
point of view of energy. Therefore there must be some factors which oppose 
the increase in entropy (i.e. folding) and encourage on the other hand the 
folding of the molecules. Since covalent structure is conserved during unfold- 
ing, forces which oppose the unfolding of the protein could not be covalent 
interactions. Transformation of the non polar molecules into water becomes 
with a considerable decrease in entropy. The entropy of hydration of those 
non polar groups must be temperature dependent since hydration leads to 
an important increase in the heat capacity. The same thing is expected for 
the enthalpy of hydration. In fact, if we look at the not hydration effect 
during the transition of the non polar components from gas phase to water 
we observe that both enthalpy and entropy are negative and the magnitudes 
decrease with increasing temperature. Hance, one could consider that the hy- 
dration effect is responsible from the dissolution of the non polar compounds 
in water. 
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